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Abstract 

We study the relations between massless Dirac fermions in an electromagnetic field and atoms 
in quantum optics. After getting the solutions of the energy spectrum, we show that it is possible 
to reproduce the 2D Dirac Hamiltonian, with all its quantum relativistic effects, in a controllable 
system as a single trapped ion through the Jaynes-Cummings and anti-Jaynes-Cummings models. 
Also we show that under certain conditions the evolution of the Dirac Hamiltonian provides us 
with Rashba spin-orbit and linear Drcsselhaus couplings. Considering the multimodc multiphoton 
Jaynes-Cummings model interacting with TV modes of electromagnetic field prepared in general 
pure quantum states, we analyze the Rabi oscillation. Evaluating time evolution of the Dirac 
position operator, we determine the Zitterbewegung frequency and the corresponding oscillating 
term as function of the electromagnetic field. 
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1 Introduction 



In graphene [H [2j [3] , the linear electronic band dispersion near the Dirac points gave rise to charge 
carriers (electrons or holes) that propagate as if they were massless fermions with speeds of the order 
of 10 6 m/s rather than the speed of light 3 x 10 8 m/s. Hence, charge carriers in this structure should 
be described by the massless Dirac equation rather than the usual Schrodinger equation. The physics 
of relativistic electrons is thus now experimentally accessible in graphene based solid-state devices, 
whose behavior differs drastically from that of similar devices fabricated with usual semiconductors. 
Consequently, new unexpected phenomena have been observed while other phenomena that were well- 
understood in common semiconductors, such as the quantum Hall effect [HE]. Thus, graphene devices 
enabled the study of relativistic dynamics in controllable nano-electronic circuits (relativistic electrons 
on-a-chip) and their behavior probes our most basic understanding of electronic processes in solids. 
It also allowed to establish a relationship between massless Dirac fermions in two dimensions and 
quantum optics, for instance one can see [DJ E]. 

On the other hand, Jaynes-Cummings (JC) model is describing the basic interaction of a two-level 
atom and quantized field, which is also the cornerstone for the treatment of the interaction between 
light and matter in quantum optics [8]. It can be used to explain many quantum phenomena, such 
as the collapses and revivals of the atomic population inversions, squeezing of the quantized field 
and the atom-cavity entanglement. Recent experiments showed that JC model can be implicated in 
quantum-state engineering and quantum information processing, e.g. generation of Fock states [9J 
and entangled states [ID] , and the implementations of quantum logic gates [IT], etc. Originally, JC 
model is physically implemented with a cavity quantum electrodynamics system, see for instance |12j . 
Certainly, there has been also interest to realize JC model with other physical systems. A typical 
system is a cold ion trapped in a Paul trap and driven by classical laser beams [13], [H] where the 
interaction between two selected internal electronic levels and the external vibrational mode of the 
ion can be induced. 

The relationship between the relativistic systems and atoms in quantum optics is studied for 
differen purposes. With this respect, the first theoretical proposal was done by Lamata et al. |15j 
and its experimental realization was done by Gerritsma et al. [16j . Subsequently, different progresses 
appreared on the subject [T7] HS1 H21 EDJ where for example in [2D], the dynamics of the 2+1 Dirac 
oscillator exactly was studied and spin oscillations due to a Zitterbewegung of purely relativistic origin 
was found. An exact mapping of this quantum-relativistic system onto a JC model is established. This 
equivalence allowed to map a series of quantum optical phenomena onto the relativistic oscillator and 
vice versa. A realistic experimental proposal was made, in reach with current technology, for studying 
the equivalence of both models using a single trapped ion. 

Motivated by different developments and in particular [15} [20] , we consider massless Dirac fermions 
in an electromagnetic field and show its link with JC and anti-Jaynes-Cummings (AJC) models. This 
can be done by determining the solutions of the energy spectrum and mapping the corresponding 
eigenspinors in terms of two states of JC model. Furthermore, we analyze the Rabi oscillations by 
considering atoms coupled to the field states and evolve in time according to the stationary eigen- 
spinors. This allows us to obtain two normalization coefficients in terms of the Rabi frequency, which 
verify the normalization condition. By evaluating the probabilities of the atom in the ground and ex- 
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cited states, we show that the probability predicts sinusoidal Rabi oscillations, much like the classical 
case, and exists even when the field is initially in the ground state. 

Finally, we deal with another aspect of Dirac fermions in graphene [2H |22| [23] that is the Zitter- 
bewegung (ZB) effect. Indeed, by calculating time evolution of the Dirac position operator, we end up 
with an extra term. After inspection, we show that this term represents a quantum oscillation about 
the classical motion and indicates the existence of the ZB effect. By evaluating the expectation value 
of the Dirac position, we obtain the ZB frequency in terms of the energy levels and magnetic field, 
which disappears for a null energy. 

The present paper is organized as follows. In section 2, we consider the Hamiltonian describing one 
massless Dirac fermion in electromagnetic field with Zeeman effect. By choosing the Landau gauge, 
we obtain the eigenvalues and the corresponding eigenspinors. After setting different ingredients, we 
establish the link between our system and quantum optics through JC and AJC models in section 
3. The Rabi oscillations for atoms coupled to field states in graphene will be analyzed in section 4. 
Using the Heisenberg picture dynamics we study the ZB effect and determine its oscillation frequency 
in section 5. Finally we close our work by concluding and giving some perspectives. 



2 Dirac fermions in electromagnetic field 

In graphene, the two Fermi points, each with a two-fold band degeneracy, can be described by a low- 
energy continuum approximation with a four-component envelope wavefunction whose components 
are labeled by a Fermi-point pseudospin = ±1 and a sublattice forming an honeycomb. Under this 
approximation we can describe our charge carriers, in electromagnetic field (E,B) where E = Ee y , 
by a gauge invariant Dirac Hamiltonian in 2+1 space-time dimensions with minimal coupling to the 
vector potential and Zeeman effect as follows 

H D = v F a ■ tt + a y eEy + ^gfi B Ba z (1) 

where vf ~ 10 6 m/s is the Fermi velocity, the conjugate momentum is tt = p — -A, g is the Lande 
factor, [ib is the Bohr magneton, a are Pauli spin matrices and the nonzero matrix elements of 
a z are (t| & z |t) = +1 and (4_| cr z \\) = —1. Choosing the vector potential in the Landau gauge 
A(x, y) = {-By, 0) to write (P) as 

H D = v F \ eB . eE c VF )+-gfi B Ba z . (2) 

\Px + iPy -^fy + iffy o J 2 

By setting £ = 1 + and £ = 1 — i^g, we obtain the form 




Px - ip y - p- iy \ 1 
H D = v F \ h - B + -gfi B Ba z (3) 

\ Px + *Py - jr tV J 2 

where Ib = \l |§ is the magnetic length. 

To diagonalize the above Hamiltonian, it is convenient to introduce the annihilation and creation 
operators. They are 

^= (iPx -Py- ijjft) , « f = ^= ("ift ~Py + ijrty ) (4) 
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which satisfy the commutation relation 
These can be used to write Hp as 

H D ~- 



a} , a 



I. 



(5) 



(6) 



—ifwo c a —^gfisB I 

where we have set the frequency oj c = \^2j^. 

The complete solutions of the energy spectrum can be obtained by solving the time independent 
Dirac equation Hd4> = £ 4>- Note that, the form of our Hamiltonian suggests to separate variables by 
writing the eigenspinors as 



{x,y) = x{x)ip{y), tp{y) = Oi,^)* , x(x) = e 



1 A." x <.£ 



(7) 



with k x is a real parameter that stands for the wave number of the excitations along x-direction. Thus, 
the reduced eigenvalue equation is given by 



fa) 



which leads to the relations 



1>i) 

fa) 



-ihio c a | Vi) = ( e + -QHbB ) | 'fa). 



(8) 

(9) 
(10) 



As we will see soon these will be solved by making separation between the positive and negative energy 
solutions as well as discussing the relation between them. 

The positive energy solution can be obtained from ([9]) and (|10p to get the following relation between 
spinor components 

— iftwc 



fa 



-a | i/ji). 



(11) 



(e + \g\iBB) 

Clearly, it is necessary to impose the condition e ^ — ^gfisB < 0, which tells us that (jlip is valid only 



for positive energies corresponding to the positive eigenspinors denoted by | ip 
from (jlip we show that | t/)f) satisfies the second differential equation 

1 V 

-gUBBj -h 2 u 2 c a)a 



4>t 



. Now, 



V>+> = (12) 
which is nothing but a simple harmonic oscillator type. Thus, the positive eigenvalues are given by 

4 = +J^n + i( W Bf, n = 0,l,2,--- (13) 



where the corresponding eigenstates take the form 



,.tr 



0). 



(14) 
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At this level, it is interesting to mention that these eigenstates are clearly depending on the electro- 
magnetic field. To clarify this statement, we write them as 

l * + >-^r(*- + *-£ & )" |0) (I5) 

with k x is the wave vector along x-direction and the parameter £ = 1 — i-^g is a function of the ratio 
between the two fields. Now inserting (|14p into (|lip to find 



I V^~) = + , ^ C a\n) = -i\\ n + 2 ^ bB I n - 1). (16) 
Combining all, we write the eigenspinors as 




u+ | n) 
| n - 1) 



z9HbB 



(17) 



(18) 



As far as the negative energy solution is concerned, we use similar approach to end up with the 
eigenvalues 



tPu2(m + l) + l(gii B B) 2 , m = 0,l,2,---. (19) 
To derive the corresponding negative eignespinors, we use the energy conservation 

(4) 2 = {e- m ? (20) 
which tells us that the two quantum numbers are related by n = m + 1. Thus, we obtain 

V \ I") / 

where it^ are given in (|18j) . We recall that all eigenspinors obtained so far are electromagnetic field 
(E, B) dependent. 

For later use, it is convenient to write the eigenspinors (|17p and (|2ip in an appropriate form. This 
can be done by using the Pauli spinors, such as 

I ) =1 ^n,a,e+) = U+ \ U^) - iu~ \ U - 1, 1) (22) 
I ^n) =\ ^n,a,e-) = U n I n i t) + *«+ | U - 1,4.) . (23) 

This summarizes the obtained results so far, which will be employed to study different issues in the 
forthcoming anaylsis. More precisely, we will implement them to establish a link between massless 
Dirac Fermions in electromagnetic field and atoms in quantum optics. Subsequentely, we will discuss 
the possibility to make contact with the Rabi oscillations and ZB effect as well. 
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3 Link with Jaynes— Cummings models 



As claimed before, we would like to use our results in order to make contact with quantum optic 
systems through two relevant models, i.e. JC and AJC. In doing so, first let us prepare our system to 
make the required contact easy and clear. From ([9]) and (|10p we show that the Hamiltonian ([6]) takes 
the form 

H D = ihw c (a) | ■0 2 )(V'i I ~ a I ^i)("02 |) + -g^ B Ba z (24) 
which can easily be written as 

H D = h {g'o'a) + g'*a + aj + ^g^i B Ba z . (25) 

This Hamiltonian can be interpreted as a model of quantum optics where g' = iuj c is the coupling 
strength between the atom and electromagnetic field, a + and cr~ are the spin raising and lowering 
operators, respectively. More precisely, ([25]) is describing a single two-state atom, represented by the 
Pauli matrices, interacting with a (single-mode quantized) electromagnetic field. This statement will 
be clarified in the next. 

We will show how to implement the dynamics of (|25p in a single ion inside a Paul trap with two 
frequencies v x and v y . For this, let us consider a system of two-level atom of the ground | e) and 
excited | /) states where the eigenspinors (|17p can be split into 

| i;) =| ^) | e)+ | ih) | /)■ (26) 

These two metastable internal states | e) and | /) are corresponding to the following energies 



Hu ne = hyJu 2 c n e + {l/2gii B B/h) 2 , e n , = hu nf = Ju*n f + (l/2gn B B/h) 2 (27) 



where the energy difference between them is characterized by the transition frequency ujq = (e ne — e nf ) / h. 
All observable quantities of this two-state system can be conveniently represented by the Pauli oper- 
ators, such as 

o* =1 e >( e I - I /)(/ I, a + =\e)(f\, <r-=|/)(e| (28) 

where the expectation value of the operator a z is the atomic inversion, while a + and a" induce 
upward and downward transition respectively. It is sometimes convenient to work with the Hermitian 
combinations of a + and o~ , i.e. a x = a + + a~ and a y = i{a + — a~). Using these to write the 
Hamiltonian (|25p as 

hjj c ( cE \ 1 

Hd = V F {VxPx + o-yPy) + y-a x + —a y jy + -gfi B Ba z (29) 

which is showing its dependence on the electric field through the extra term ^-<7 y appearing in the 
Hamiltonian form. Clearly for E = 0, one recovers a system of massless Dirac fermions in magnetic 
field and thus the extra term can be seen as a shift to the spin a x . Certainly, this will play a crucial 
role in the forthcoming analysis and make difference with respect to the former results obtained in 
this direction 
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Having set all ingredients, let us show that the Dirac Hamiltonian (|29p can be linked with JC 
and AJC models. We start by introducing JC model, usually called red-sideband excitation, which is 
consisting of a laser field acting resonantly on two internal levels. Typically, a resonant JC coupling 
induces an excitation in the internal levels while producing a deexcitation of the motional harmonic 
oscillator and viceversa. The resonant JC Hamiltonian is given by 



H 



JC 



(30) 



where a, and af (i = x, y) are the annihilation and creation operators associated with a motional 
degree of freedom, rji := ki\J%j2Mvi is the Lamb-Dicke parameter [14] . ki is the wave number of the 
driving field, M is the ion mass, z/j are the natural trap frequencies, <j) is the red-sideband phase, flj 
and Si are the excitation coupling strengths. Using the realization @ for the annihilation and creation 
operators to map (|30|) as 

jc Vi^B 



H 



V2 



i(a + e i(t> - a-e- i(t> )p x - {a + e i(t> + a'e'^Py 

cE 



-i(a + e l4 



a e 



V2l B 

On the other hand, the AJC Hamiltonian reads as 



t> 



y + HSi(T z . 



H 



AJC 



where 99 is the blue-sideband phase. In similar way to (|30p . we show that Hf JC takes the form 



H 



AJC 



[-i{a + e iv - a-e~^) Px - (a + e itfi + a~ e^py] 



(31) 



(32) 



(33) 



+ 



V2l B 



i(a e 



up 



a e 



cE 

— (a+e^ + a-e-^) 



By fixing different parameters, let us see how does look like the sum of Hf c + Hf Jc . Indeed, for 
particular values of the involved parameters we obtain table 1: 



i 


Si 


4> 


<P 


H JC + H AJC 


X 


5 


3tt 
2 


7T 

2 


\f2rfll B u x p x + h5a z 


y 





3tt 
2 


3tt 
2 


\f2rfll B a y Py + °§ a y y 


y 





3?r 
2 


7T 

2 





Table 1: Linear combinations of JC and AJC Hamiltonian's for given values of the involved parameters 

Now by summing all terms in table 1, one gets 



{H? C + Hf JC ) = V2i]Ql B (a xPx + Oypy) + ^ ?? ° f-a x + ^y^j V + ^5a z . 



(34) 



A comparison to the Dirac Hamiltonian (|29p yields to the following correspondence 



where f2 = Jlj and r\ = rji, Vi = x,y. (j35|) means that if the Lamb-Dicke parameter r/ is restricted 
to the ratio between the cyclotron frequency co c and the excitation coupling strength supplemented 
by a variation of magnetic field like the excitation coupling strength 5, then the present system can 
be linked with the atoms in quantum optics. Otherwise, (j35[) tells us that it is possible to reproduce 
the 2D Dirac Hamiltonian (|29|) . with all its quantum relativistic effects, in a controllable system as a 
single trapped ion. This conclusion is deduced also by analyzing massive Dirac fermions in magnetic 
field in the same framework |20j. 

We close this section by showing that the system under consideration can be reduced to spintronic 
ones. In doing so, let us choose different parameters as summarized in table 2: 



i 


6 


4> 




H JC + H AJC 


X 


5 





7T 


\f2rfll B o y p x + Ma z 


y 











-^2rini B a xPy -^Ml§a x y 



Table 2: Configuration of four parameters leading to new form of H( + . 

Summing over i to obtain the Hamiltonian 

V (Hf c + Et JC ) = V2 V m B (a yPx - a xPy ) - ^^ C -§a x y + H5a z . (36) 

By taking E = and 5 = in ()36p . the resulted term is nothing but the Rashba spin-orbit coupling [23] 
where the coefficient ^/2?]QIb is identified to the coupling parameter of Rashba type. 

One can ask for other configurations of the involved parameters and see what we gain. Indeed, let 
us take the choice listed in table 3: 



i 


5 





f 


H JC + H AJC 


X 


5 


3tt 
2 


7T 

2 


V2r)m B (TxPx + frScfz 


y 





7T 

2 


7T 

2 


\f2r]VLl B (TyPy c B o y y 



Table 3: Configuration of four parameters leading to new form of H( + Hf . 

As before we sum up to obtain the Hamiltonian 

£ (H t JC + Ht J °) = V^flifl (cr xPx - oypy) - ^121^ + Ua z (37) 

i 

which reduces to the linear Dresselhaus coupling [23] for E = and 5 = 0. According to ()36p and ([3711 . 
one can conclude that the system under consideration is also sharing some features with spintronic 
systems. Therefore one can go further to investigate the basic features of (|36|) and (|37|) to extract 
more information about these relationships. 



4 Rabi oscillations 

The Rabi oscillations of fermions in graphene have been studied in different contexts to deal with some 
issues [Q[7]. For example in [U] the relation between the canonical model of quantum optics, JC model 
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and Dirac fermions in quantizing magnetic field was investigated. It was demonstrated that Rabi 
oscillations are observable in the optical response of graphene, providing with a transparent picture 
about the structure of optical transitions. While the longitudinal conductivity reveals chaotic Rabi 
oscillations, the Hall component measures coherent ones. This opens up the possibility of investigating 
a microscopic model of a few quantum objects in a macroscopic experiment with tunable parameters. 

Motivated by the above developments, we analyze the Rabi oscillations of the system under con- 
sideration. We start by expanding a general state on to the obtained eigenspinors ()17p to study the 
corresponding dynamics. Let us assume that initially at t = the atom is in the lower state | /) 
and the cavity contains precisely n — 1 photons, i.e. the field is in a number (Fock) state | n — 1) 
with n = 1, 2, • • • . Then the initial state of the system is | f,n— 1), while the interaction term of the 
Hamiltonian (|25p connects this initial state to the final state | e, n). To proceed further, we consider 
the initial state for the field as | ip(0)) =| f,n — and assume an atom in the excited state is 
injected into the field. Using (|22|) . (|23p and (|26|) . one can write the initial state as 

l^(0)>=i«" |^ w +>-i«+|^ w -). (38) 

Its evolution state can easily be obtained by acting as 

\^{t))=e-r^ t | V(0)) =C e , n (t) | e,n,t) + C /in _i(t) | f,n-l,l) (39) 

where the coefficients are given by 



/ 4Cn 

Ce,n (t) = J Y^^l Sill^nt) (40) 

i 

C/ „_i(t) = cos(u) n t) + 7== sin(w n i). (41) 

VI + %n 

It is easy to check that they verify the normalization condition \C e ^ n (t)\ 2 + |Cj jn _i(£)| 2 = 1. The 
frequency of oscillations is given by 

w„:=y = <Vl + 4Cn (42) 

where the parameter Q = (|^) 2 controls the nonrelativistic limit. (|39p shows very clearly the non- 
stationary states of localization (| e, n, f), | f,n — 1, \)), which is nothing but the Rabi oscillation that 
is a simple quantum beat interstate | e, n, "f) and | f,n — l,\). Moreover, the dynamics of (|39p is 
completely similar to the atomic Rabi oscillations that can be seen from the JC model. This similarity 
shows another way how the system under consideration can be linked with Rabi oscillations. 

To understand better the present situation, let us analyze the probabilities of existence. Indeed, 
for a atom in the state | f,n— 1,4-), we find 

f , s 4£n 9 , , / S 2 \ 9 , 



PLlit) = 1 " T^i^ = 1-11-^) sin^a;„t) 



where uj n is now identified to the Rabi frequency cjr. It is clear that, -P r f_ 1 (t) oscillates between a 



maximum of unity and a minimum of ( — ) . ()43|) predicts sinusoidal Rabi oscillations much like the 
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classical case. However, even when the field is initially in vacuum state, i.e. n = 0, the probability 
exists. 

Now let see what happens when the atom is in the state | e,n,\). Indeed, after some algebra we 
get the probability 

Pn(t) = sin 2 l (uRt) = (l " |r) smW). (44) 

This can be plotted as shown in Figure 1 for some particular value of 5. One can easily notices that 
for the resonance case (5 = 0), P%(t) reduces to 

P*(t) = sm 2 (oj R t) (45) 

which is in agreement with blue figure (Figure 1). For t = tt/2ur all the atomic population has been 
transferred to the excited state. Clearly, the probabilities P£(t) and -P/_ 1 (i) oscillate in time with the 
Rabi frequency ujr where their amplitude are maximum at 5 = and decrease rapidly as 5 increases. 
Contrary, in a stationary state the probabilities of localization of the particle are equal and constant. 




co t 

n 

Figure 1: Plots of P^(t) versus t for various detunings 5. 

It is frequently convenient to consider the atomic inversion W(t), which is defined as the difference 
in the excited and ground state populations. This is 

W(t) = P n e (t) - PU(t) = 2 (l - 4) sinW) - 1 (46) 

which can be simplified by assuming that the atom is initially in the ground state and requiring the 
resonant case. These allow us to obtain 

W(t) = sin 2 (u R t) - cos 2 (w R t) = - cos(2w R t) (47) 

where the Rabi frequency reduces now to the quantity uir = y/nuj c . Again, for t = 7r/2y / nu; c all the 
population is transferred to the excited state and we have W(n /2y/noj c ) = 1. On the other hand, if 
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we choose t = n/4^/nu; c , then W (tt / '4y/nu) c ) = and the population is shared coherently between the 
excited and ground states. This case shows the constraint 

C e , n (7r/4Vnw c ) = C /jn _i(vr/4V^w c ) = (48) 

and therefore the corresponding state is 

IV^/V^c)} = 4| (|e,n,t) + |/,n- U» • (49) 

In summary, one can see that the behavior of two-level atom interacting with a single electromagnetic 
mode is surprisingly rich. In fact, there are Rabi oscillations that collapse and remain quiescent, 
revive, then collapse again. 



5 Zitterbewegung effect 

The Zitterbewegung (ZB) effect, called also trembling motion, is consisting of a helicoidal motion of 
a free Dirac particle. It is a natural consequence of the non-commutativity of its velocity operator 
components VF&j, with i = x,y. In connection with quantum optics realm, this quantum relativistic 
effect was studied by Lamata et al. [15]. The ZB effect has been also studied for different reasons in 
systems made of graphene |2 H I22 ( |2"3"]. For instance in |22j . the electric current and spatial displacement 
due the ZB effect of electrons in graphene in the presence of an external magnetic field were described. 
It was shown that, when the electrons are prepared in the form of wave packets, the presence of a 
quantizing magnetic field B has very important effects on ZB. (1) For B ^ the ZB oscillations are 
permanent while for B = they are transient. (2) For 5^0 many ZB frequencies appear while 
for B = only one frequency is at work. (3) For B ^ both interband and intraband (cyclotron) 
frequencies contribute to ZB while for B = there are no intraband frequencies. (4) Magnetic field 
intensity changes not only the ZB frequencies but the entire character of ZB spectrum. An emission 
of electromagnetic dipole radiation by the trembling electrons was proposed and described. It was 
argued that graphene in a magnetic field is a promising system for an experimental observation of 
Zitterbewegung. 

Let us study the ZB effect for massless Dirac fermions in an electromagnetic field and discuss what 
will add an electric field to the results presented in |22j . In doing so, we describe the dynamics of the 
system under consideration through the Heisenberg picture. Indeed, we calculate the time evolution 
of the Dirac position operator 

x(t) = e~Tr HDt x e^ Hot (50) 
which satisfies the Heisenberg equation of motion 

= -*i H D, x] = vpaj (51) 

where Hp is the Hamiltonian given in (TjQ). To get the explicit form of x(t), we introduce the dynamics 
of the Pauli operators 

-ih^f = [H D ,a j ] (52) 
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which can be evaluated by using the anti-commutation relations {ai,aj} = 2<5jjl = x,y,z). This 
leads to the equation 

- ifr^jr = -2v F TTj - 2eEy5 yj + 2H D dj. (53) 
It is convenient for later use to introduce the compact form 

- ih^- = 2H D ff (54) 

where the operator rj is given by 

fj = dj — vpH^TTj — eEyHp Syj (55) 
and satisfies the following dynamics 

It can be solved easily to find 

ff(t) = e T H ^ff (57) 

where ffo is a constant operator 

Vo = V(0) = ffj (0) - v F Ho X T?j ~ eEyH D l 5 y3 . (58) 
Since we have {H£>,ff\ = = {Hp, ffo}, thus (f5T|) rewrites as 

ff(t) = Voe^ HDt . (59) 



Combining (|53j) . (j55|) and (|59j) . to end up with 

^ = vlH^TTj + eEvpyH^Syj + v F ffoe~^ HDt (60) 

which can be integrated to get the time evolution 

x(t) = vpH^TTjt + eEvpyH^Syjt + -ihv F ffoH^ x e~^ HDt + x(0) (61) 

where x(0) is a constant (operator) of integration 

^(0) ~ IjifrvFVjifyHp 1 + ^ihv F Hp 2 irj + -ihv F eEyH^ 2 5 y j. (62) 

Now using the correspondence (|35p . which is linking our system with atoms in quantum optics, we 
map (|6"Tj) as 

m = m + 2 V %n% + V2eErjl B ny Syjt 

( _ . y/2rjl B ^j + eEy5 vj \ i\/2hrjl B VL / 2*^ t \ 

+ h (t) ir D — -)-2H^-\}- eh )■ 
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This result is showing different contributions originated from different sources. Indeed, the first two 
terms on the r.h.s. account for the classical kinematics of a free particle in electromagnetic field (E, B). 
However, the last term represents a quantum oscillation about the classical motion with a frequency 

1\E I 

1 h Dl where Ed = e n . This oscillating term is the so-called ZB effect or "trembling motion" and its 
frequency associated with the measurable quantity (x(t)) can be estimated as 

u ZB = Hl^i = 2 v / VQ 2 |n| + <5 2 (64) 

where the correspondence (|35p is used. On the other hand, the presence of the electric field E is 
clearly shown in equation ([63]) . Thus, one can interpret the terms (V2eEr]lB^lyS y j) and (eEy5 y j) as 
shift with respect to the original terms, which of course disappear if we require E = 0. This tells us 
that our results are interesting and generals such that one can recover the standard ones by simply 
fixing some parameters. 



6 Conclusion 

The solutions of the energy spectrum of a system of massless Dirac fermions in the electromagnetic 
field are obtained. The corresponding Dirac Hamiltonian Hp is used to investigate two-level atom 
interacting with a single electromagnetic field. The connection between JC model in quantum optics 
with Hp is established in a simple manner without taking any limit on the strength of the electro- 
magnetic field. This can be served as a tool to study atomic transitions in quantum optics using the 
relativistic quantum mechanical models in the presence of the electromagnetic field. 

More precisely, after mapping the eigenspinors in terms of ground and excite states of two-level 
atom, we derived a new form of the Dirac Hamiltonian. This is used to establish a link with JC and 
AJC models as linear combination under the correspondence between their parameters. Making a 
suitable choose of the involved parameters in the game, we showed that the linear combination can be 
reduced either to the Rashba or Dresselhous couplings. This is a significant results because it relies 
two sectors such graphene and quantum optic systems. 

Because of the Rabi oscillations are occurring in the JC model and its relation to our system, we 
analyzed them. In fact, we mapped the eigenspinors in terms of two level atom states, i.e. ground and 
excited and showed that the Rabi frequency can be obtained in terms of the quantized energy levels. 
Calculating the probabilities for atom in ground and excited states, separately, we obtained different 
results showing that these probabilities are maximum when the detuning 5 is null and decrease rapidly 
as 6 increases. These are used to discuss the resonant case by evaluating the atomic inversion W(t) 
and conclude that the population of the state is depending on the value attributed to W(t). 

Finally, we analyzed the Zeitterbewegung effect for massless Dirac fermions in the electromagnetic 
field. This was done by using the Heisenberg picture dynamics to evaluate the time evolution of the 
Dirac position operator. This allowed us to obtain a quantum oscillation about the classical motion 
with a defined frequency. Its oscillation between the negative and positive energy states is similar to 
the Rabi oscillation in the two-level atom. By evaluating the position average, we obtained the ZB 
frequency in terms of the quantized energy system. 

We close this part by listing some ideas, which can be worked by adopting our approach presented 
in this paper. Indeed, it will be of interest to analyze the massless Dirac fermions in the presence of 
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a perpendicular magnetic field and by considering an electric field supported by 2 x 2 unit matrix as 
formulated for graphene in [3]. This certainly will bring new results and interesting conclusions due 

to the reduction of the magnetic field from B to B' = B^ 1 — (^y^j^J ■ One also can study in this 
framework the tilted magnetic field case instead of an uniform one. 
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